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Abstract 

The time evolution of a Gamma-ray burst (GRB) is associated with 
the evolution of a supernova remnant (SNR). The time evolution of 
the flux of a GRB is modeled introducing a law for the density of the 
medium in the advancing layer. The adopted radiative model for the 
GRB in the various e.m. frequencies is synchrotron emission. The X-ray 
ring which appears a few hours after a GRB is simulated. 
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1 Introduction 

Gamma-ray bursts (GRBs) started to be observed from an astronomical point 
of view by [TJ. Recently some excellent reviews have been written, see [21 
El IH IS]- The recent trend of research connects GRBs with supernovae (SN) 
explosions and their consequent remnants (SNRs). The afterglow is observed 
in X, optical, and infrared, see [B]. The established connection between GRBs 
and SNRs allows of restricting the astrophysical environment to the advancing 
shock. In the last 60 years the theoretical studies of the advancing shock of 
SNR have been focused on the relationship with the instantaneous radius of 
expansion, R, which is of type oc t m , where t is time and m is a parameter 
that depends on the chosen model. The most popular model is the Sedov- 
Taylor expansion which predicts R oc t 0A , see [H EJ [9j HU] . and the thin layer 
approximation in the presence of a constant density medium, which predicts 
R oc t 25 , see [11]. A simple version of the Sedov-Taylor solution, which also 
considers the transition from the relativistic to the non-relativistic regime, 
appeared in [T21 [2]. The radiative transfer during the afterglow phase has 
been analysed by [15] . where a comprehensive analysis, including images of the 
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expected light curves, is carried out. The actual situation of the theoretical 
and observational research on SNRs leaves a series of questions unanswered or 
only partially answered: 

• Can we deduce a theoretical expression for the counts/s versus time in a 
GRB assuming a direct proportionality between the flux of kinetic energy 
of the advancing shell and the intensity of emission? 

• Can we simulate the X-ray ring of GRB 032103, which was observed six 
hours after the burst? 

• Can we deduce a theoretical expression for the frequency at which the 
GRB becomes visible function of the time? 

In order to answer these questions, Section [2] reports the data on GRBs. Sec- 
tion E] reviews the free expansion, the Sedov-Taylor solution and the power law 
model. Section H] reviews the existing situation about the radiative transport 
equation and a simple geometrical model which leads to the ring emission. 
Section [5] contains some scaling arguments that lead to the observed count/s 
versus time relationship. 

2 The Observations of GRBs 

GRBs are characterized by a rapid rise (called the trigger) in the rate of de- 
tected photons per second at a given energy, for example 20 kev. A typical 
example observed by the Burst and Transient Source Experiment (BATSE), 
see [14] . on board the Compton Gamma Ray Observatory (CGRO) is reported 
in Fig. [TJ 

^From the previous figure it is important to point out that the time nec- 
essary to reach the maximum is 1.25 s. In the following we will consider this 
time small compared to the time in which the counts/s decreases. 

A numerical analysis of the observed rate versus time relationship can be 
done by assuming a power law dependence for the rate of counts of the type 

rate(t) = Ct 5 counts s' 1 , (1) 

where the two parameters C and 5 are found from a numerical analysis of the 
data. 

An interesting analysis of the time evolution of 12 GRBs as observed by 
the Swift Gamma- Ray Burst Mission (Swift), see [15J, is reported in [16J. We 
selected one GRB from the previous list, namely, GRB050814 XRT 1.5-10 kev 
which is visible in Fig. 3 (top-right) of [16J. For the previous GRB as well the 
fit with relationship (JTJ) is reported in Fig. [2j 
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Figure 1: Light curve of BATSE trigger 8121 in the energy channels 1- 
4 ( > 20 kev) (full line) and power law fit (empty stars). The real data 
have been extracted by the author from the BATSE database available at 
the web site http:/ /www.batse.msfc.nasa.gov/batse/grb/light cilTve/^[ In this 
case 5 = -0.1086 and C = 20382.27 when t > 1.25 s. 
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Figure 2: Light curve of GRB050814 XRT 1.5-10 kev (empty stars) and 
power law fit (dotted line). The real data have been extracted by the author 
from Fig. 3 (top-right) of Willingale et al. 2010. In this case 5 = —1.22 and 
C = 4214. 
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Figure 3: Count spectrum in energy as measured by the three instruments 
on board the Fermi Gamma-ray Space Telescope, which are Nal, BGO, and 
LAT. The data are extracted by the author from Fig. 2A in [17]. 

The count spectrum in energy of GRB 080916C, see [TF], is reported in 
Fig. [3] and can be approximated by the following formula 

E(kev) « 1563E(kevy 1 ' 7 when 10 kev <E< 10 7 kev. (2) 

Another interesting property of the GRBs is that the frequency of observa- 
tion at which they start to be visible decreases with time. An example of such 
behaviour for GRB 050904 in the BAT, XRT, J, and I bands can be found in 
Table 1 of [18] . The numerical analysis of the data in Fig. H] gives 

v = 1.0310 25 t(s)- 2 - 27 Hz (3) 

where t has been chosen as the value of time at which the flux at the chosen 
frequency is maximal. 

Also interesting is the case of GRB 032103, which was observed in the X-ray 
(energy range 0.7-2.5 keV) by XMM-Newton's MOS cameras six hours after 
the GRB, see [19]. In this case the combined MOS images show two concentric 
rings which expand outwards. The numerical analysis gives 9 oc (t — to) 1 / 2 
where 9 is the angular separation between the the two maxima in the radial 
profiles of counts, t is the time, and t the initial time. In the following we will 
associate the GRBs with the explosions of supernovae which are classified as 
SNR, see Section 4 in @]. 



running head 




Log (t [sec]) 



Figure 4: Bands from the radio through the IR/optical to X-ray and BAT at 
which the expansion becomes optically thin. The energy bands are expressed 
in logarithms of the frequency in Hz along the y-axis. The x-axis reports the 
logarithm of the time in seconds. The data are extracted by the author from 
Table 1 in [IB]. 
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3 Some existing solutions 

This section reviews the expansion at constant velocity, the Sedov-Taylor, and 
the power law models. 



3.1 The constant expansion velocity model 

The SNR expands at a constant velocity until the surrounding mass is of the 
order of the solar mass. This time, t^, is 

t M = 186.45 1. ° yr , 4 
ynov 100 oo 

where M & is the number of solar masses in the volume occupied by the SNR, n , 
the number density expressed in particles cm -3 , and Vioooo the initial velocity 
expressed in units of 10 000 km/s, see [20j. 

3.2 The Sedov-Taylor solution 

[3] The Sedov-Taylor solution is 

^1/5 
4 7T p 



„. . /25 Et 2 . 

R (t) = hr — ■ ( 5 ) 



where E is the energy injected into the process and t is time, see [8j |9j [T0~| 120]. 
Our astrophysical units are: time, (ti), which is expressed in years; E$\, the 
energy in 10 51 erg; no, the number density expressed in particles cm -3 (den- 
sity p = n m, where m = 1.4m H ). In these units, Equation (J5J) becomes 



R(t)n 0.313 pc . ( 6 ) 
V n 

The Sedov-Taylor solution scales as t 0A . We are now ready to couple the 
Sedov-Taylor phase with the free expansion phase 

f 0.0157 tpc ift<2.5yr 
{> \ 0.0273 v^pc ift>2.5yr. {> 

This two-phase solution is obtained with the following parameters M Q = 1, 
n = 1.127 10 5 , E51 = 0.567 and Fig. [5]reports its temporal behaviour as well 
as the data. ^From the previous figure, it is clear that the free expansion plus 
the Sedov-Taylor phase is not a satisfactory model. 
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Figure 5: Theoretical radius as given by the two-phase solution (full line), 
and astronomical data of SN 1993Jwith vertical error bars. 

3.3 The power law model 

We shall discuss a power law dependence of the type 

R(t) = Ct a , (8) 

where the two parameters C and a can be found from the observations. Ten 
years of observations of SN 1993 J , see [21], allow of fixing these two parameters: 
a = 0.828 ± 0.0048 and C = 0.0155 ± 0.00011, and Figure © displays the data 
as well as the fit. 

This observed relationship allows of expressing, from an astrophysical point 
of view, the radius as 



R(t) = 0.0155 1 



0.828 



pc 



(9) 



where the time t is expressed in years. The velocity in this model is 

V{t) = Ca&- 1) (10) 

and the astrophysical version 

km 



V(t) = 12587.67t 



-0.171 ' 
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Figure [7] reports the observed instantaneous velocity as deduced from the finite 
difference method and the best fit according to (jTTj) . 
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Figure 6: Theoretical radius as given by the power law model (full line), and 
astronomical data of SN 1993J with vertical error bars. 




Figure 7: Theoretical velocity as given by Eq. (ITTj) (full line) and instanta- 
neous velocity of SN 1993J with uncertainty. 
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4 Radiative processes 

This section reviews the existing knowledge about the transfer equation and a 
simple model which produces an X-ring. 



4.1 The transfer equation 

The transfer equation in the presence of emission only, see for example [22] or 

[23], is 

= -k u pl v + j v p , (12) 

where I v is the specific intensity, s is the line of sight, j v is the emission 
coefficient, k u is a mass absorption coefficient, p is the density of mass at 
position s, and the index v denotes the frequency of emission of interest. The 
intensity of radiation, i.e., the photon flux, is here identified with the counts 
at a given energy. The solution to Equation (fl2"|) is 

lM = ^(l-e~^) , (13) 

where r v is the optical depth at frequency v. 

dr v = k u pds . (14) 



The framework of synchrotron emission, as described in sec. 4 of [24J, is often 
used in order to explain a GRB, see for example [23 [251 [23 EE]- The volume 
emissivity (power per unit frequency interval per unit volume per unit solid 
angle) of the ultrarelativistic radiation from a group of electrons, according to 
EH, is 



e(u) = I P(v)N(E)dE , (15) 

where P{v) is the total power radiated per unit frequency interval by one 
electron and N(E)dE is the number of electrons per unit volume, per unit 
solid angle along the line of sight, which are moving in the direction of the 
observer and whose energies lie in the range E to E + dE. In the case of a 
power law spectrum, 

N(E)dE = KE^dE , (16) 

where K is a constant. The value of the constant K can be found by assuming 
that the probability density function (PDF, in the following) for relativistic 
energy is of Pareto type as defined in [30J : 

f{x;a,c) =ca c x~ (c+1) , (17) 
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with c > 0. In our case, c = 7 — 1 and a = E m i n , where E m i n is the minimum 
energy. /,From the previous formula, we can extract 



k = iv ( 7 - i)e; 



7-1 

vain 1 



where N is the total number of relativistic electrons per unit volume, here as- 
sumed to be approximately equal to the matter number density. The previous 
formula can also be expressed as 

*=t/— (7-l«£ , (19) 
1.4 mn 

where m# is the mass of hydrogen. The emissivity of the ultrarelativistic syn- 
chrotron radiation from a homogeneous and isotropic distribution of electrons 
whose N(E) is given by Equation (jl~6l) is, according to [29] . 



3vP = (20) 

0.933 x 10- 2 V(7)™l 7+1)/2 ( 6 - 26Xl ° 18 ) ( ^ 1)/2 

v 

ergsec -1 cm _3 Hz -1 rad -2 , 



where v is the frequency and oll{i) is a slowly varying function of 7 which is 
of the order of unity and is given by 

*M = 2 (- 3 >/ 2 l±^r(^)r(^) , (21) 

l\d 7 + 1 v 12 ; v 12 ; ' K J 

for 7 > \. 

We now continue to analyse the case of an optically thin layer in which r v 
is very small (or k v is very small) and where the density p is replaced by the 
concentration C(s) of relativistic electrons: 

j v p = K e C{s) , (22) 

where K e is a constant function of the energy power law index, magnetic field, 
and frequency of e.m. emission. The intensity is now 

I v {s) = K e C{si)dsi optically thin layer. (23) 

J SO 



The increase in brightness is proportional to the concentration integrated along 
the line of sight. 
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Figure 8: The two circles (sections of spheres) which include the region with 
constant density are represented through a full line. The observer is situated 
along the x direction, and three lines of sight are indicated. 
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4.2 The X-ring 

We assume that the number density of ultrarelativistic electrons C is constant 
and in particular rises from at r = a to a maximum value C m , remains 
constant up to r = b, and then falls again to 0. This geometrical description 
is reported in Fig. [Sj The length of sight, when the observer is situated at 
the infinity of the x-axis, is the locus parallel to the x-axis which crosses the 

position y in a Cartesian x y plane and terminates at the external circle of 

radius b. The locus length is 



l 0a = 2x^b 2 -y 2 -^/a 2 -y 2 ) ;0<y<a 

l ab = 2x ^b 2 -y 2 ) ;a<y<b . (24) 

When the number density of ultrarelativistic electrons C m is constant between 
two spheres of radius a and b the intensity of radiation is 

I 0a = C m x2x {^Jb 2 -y 2 - ^a 2 -y 2 ) ;0<y<a 

I ab = C m x2x {^jb 2 - y 2 ) ;a<y<b . (25) 

The ratio between the theoretical intensity at the maximum (y = a) and at 
the minimum (y = 0) is given by 



7<^o) = ^r • (26) 

The parameter b is identified with the external radius which denotes the tran- 
sition from the perturbed medium. The parameter a can be found from the 
following formula: 

< I(y=a) x 2 _ 
•l(y=0) >obs 



((i(J=0))L + i) 



(27) 



where ( jfe^ )obs is the observed ratio between the maximum intensity at the 
rim and the intensity at the centre. A cut in the theoretical intensity of GRB 
031203 XMM-Newton is reported in Fig. Oand a theoretical image in Fig. [TDJ 

A comparison should be made with Fig. 7 in 



5 Flux density versus time 

The source of synchrotron luminosity is assumed here to be the flux of kinetic 
energy, L m : 

L m = ^p L rcR 2 V 3 , (28) 
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Figure 9: Cut of the mathematical intensity I of the ring model, Equa- 
tion (J25D, crossing the centre (full line) of GRB 031203 XMM- Newton obser- 
vation. The x- and y-axes are in arcsec, a = 17.6 arcsec, b = 22 arcsec and 
jvyZo) = 3 as in Fig. 2 panel 6 of [T9] . 
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Figure 10: Contour map of / adjusted to simulate GRB 031203 XMM-Newton 
observation. The x and y axes are in arcsec, other parameters as in Fig. |9j 
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where R is the instantaneous radius of the SNR, and pi is the averaged density 
in the advancing layer. The averaged density in the advancing layer can be 
parametrized as 

(p^R 3 ) 1 ^ , x 

where AR is the thickness of the advancing radius in which the mass resides. 
In order to quantify the thickness we assume that 

A fl = ^i , (30) 

n 

with t]j derived from the fit with the observations and n a sufficiently large 
number such that 

/?(/ """ ) " J <i2(tm M ) , (31) 



n 

where t max is the maximum time considered. 

Therefore, the mechanical luminosity transforms to 

L m cxp- L R 3 V 3 . (32) 

The temporal and velocity evolutions can be given by the power law depen- 
dencies of Equations (JSJ) and ({TO"]) and therefore 

L m oc t~~s+p~ . (33) 

The synchrotron luminosity L\ and the observed flux S\ at a given wavelength 
A are assumed to be proportional to the mechanical luminosity, and therefore 

t 6 + 7)jp 

Sx = S (— )-— , (34) 

where So is the observed flux at t — tso- As an example from the astronomical 
data of Fig. [21 we deduce that rjj = 1.095. The progressive transparency of 
the GRB to the low frequencies with time can be explained by the dispersion 
relation of an e.m. wave in a cold plasma, which is 

uj 2 = oo 2 pe + k 2 c 2 , (35) 

where u is the frequency of the e.m. wave, u pe is the plasma frequency of the 
electrons, k is the wave-vector of the e.m. wave, and c is the velocity of light, 
see [31J . The plasma frequency of the electrons is 

u pe = (4nn e e 2 /m e ) 1/2 , (36) 



or 



f P e = w pe /27T = 8.98 x 10 3 n £ 1/2 Hz , (37) 
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where f pe is the electron plasma frequency expressed in Hz. A visual inspection 
of the dispersion relation ( 1351) allows us to say that in an astrophysical envi- 
ronment in which the number density of the electrons is decreasing with time, 
the frequency at which the medium becomes transparent will also progres- 
sively decrease. The time dependence of the frequency at which the medium 
is transparent, according to (I2"9~j) . is 

uocr 1/2 — — P — , (38) 

where T]n is the same parameter as in ( 130]) and the index II means a second 
way to deduce the free parameter. A comparison of the previous formula with 
the observational formula ([3]) gives rjn = 3.42. 



6 Conclusions 

The first two parts of the law of motion for a SNR are thought to be a free 
expansion in which R cx t and an energy conserving phase, the so called Sedov- 
Taylor phase, in which R oc t 2 / 5 . A careful analysis of SN 1993J in the first ten 
years suggests instead that R oc t 82 . 

The scaling laws of the flux of GRBs with time can be consequently de- 
duced assuming that the intensity of emission is proportional to the flux of 
the kinetic energy of the advancing shell, see Equation (1341) . This result is 
reached assuming that the thickness of the advancing layer increases oc R v . 
The comparison between the observed flux and the time of GRB050814 XRT 
1.5-10 kev gives rjj = 1.095 and the comparison between the e.m. frequency 
of transparency and time gives rjjj = 3.42. The X-ray ring connected with 
GRB 031203 XMM-Newton is simulated introducing two hypotheses: (i) the 
emission of radiation is localized in a thin layer of the advancing expansion, 
(ii) the emitting medium is supposed to be optically thin. 
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